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1. Introduction 

The present paper can be viewed as a continuation of the paper [1], where the transition 
probabihties for the totaUy asymmetric simple exclusion process (TASEP) with parallel 
update have been obtained for the ID infinite lattice, generalizing well-known result of 
Schiitz [2] for the continuous time TASEP. The analytic method, developed in the first 
half of that paper, was based on the use of Bethe eigenvectors of the evolution operator 
obtained in [3]. The problem was to find a continuous spectrum, i.e. the integration 
measure, which would allow one to construct the solutions of the master equation out of 
the Bethe eigenvectors, and in particular would give a resolution of the identity operator 
in the integral form. The proof of the formula for the resolution of the identity operator, 
was the main technical result of the paper which yielded the final determinant formula 
for the transition probability. 

A peculiarity of the finite ring is that the spectrum becomes discrete due to 
the periodic boundary conditions, being defined by the system of the algebraic Bethe 
equations (BE). These equations usually cannot be solved exactly. Most attempts to 
extract any information from them are related to the thermodynamic limit, where they 
can be reinterpreted in terms of a single integral equation. For the finite lattices, not 
only the exact form of the spectrum but even the issue of its completeness are far from 
being well understood for most of the integrablc models. 

Fortunately, for the TASEP the situation is a bit better due to a very special 
factorization property possessed by the BE. This property was observed already in one 
of the first works on the continuous time TASEP by Gwa and Spohn [4] , where it was 
used to get the asymptotical behaviour of the spectral gap. One of the most impressive 
results exploiting the special structure of the BE in TASEP was the exact derivation 
of the largest eigenvalue of the equation for the generating function of the integrated 
particle current, due to Dcrrida and Lcbowitz [5]. Using the Cauchy residue theorem, 
they managed to evaluate explicitly the sum over roots of the BE for the particular 
solution corresponding to the groundstate of the evolution operator. Recently, some 
peculiar details of the spectrum structure for the continuous time TASEP have been 
also studied in [6]. 

The solution of the Cauchy problem, i.e. finding the solution, given the initial 
conditions, for the master equation for the TASEP on the ring in continuous time, and 
also with backward ordered update, has been recently proposed by one of the authors 
[7]. To our knowledge, this is still the only example of full solution to the Cauchy 
problem for an integrable model that cannot be reduced to free fermions. The method 
of the solution is based on a geometric approach to the Bethe ansatz (BA), which treats 
trajectories of the interacting particles as free, noninteracting, but supplied with the 
additional statistical weights in such a way that the interaction is taken into account. 
The idea of the solution on a ring was to represent it as an infinite line with periodic 
patterns of synchronously moving interacting particles. As in the geometric formulation 
one works only with the ensembles of particle trajectories of a finite length, the issues 
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related to the structure of the eigenspace of the evolution operator can be by-passed. 
Motivated by the geometric solution, its analytic counterpart called the detailed BA 
has been also proposed [8]. Being defined on the infinite lattice, the detailed BA has 
a form of the infinite formal sum, and can be considered as a generating function of 
particle trajectories. As a result, the problem becomes infinite lattice-like and one can 
again use integration over the continuum to reconstruct the transition probabilities. The 
mathematical meaning of the term-by-term integration of the infinite formal series was 
not yet well understood, and, as we will see below, it is just equivalent to summation 
over the discrete spectrum given by the BE. 

The geometric method was applied to the TASEP with parallel update on the 
infinite lattice in the second half of our first paper [1]. Its generalization to the ring 
geometry follows in line with the continuous time version and as such is straightforward. 
Therefore, the aim of this paper is not only to obtain the solution of the Cauchy problem 
for the TASEP with parallel update, but also to establish a bridge between existing 
solutions on a ring and the standard BA techniques. Namely, we show that the idea of 
Derrida and Lebowitz exploiting the Cauchy residue theorem to obtain the roots of the 
BE corresponding to the groundstate solution can be developed much further and used 
to perform an exact summation over the whole spectrum of solutions. In such a way 
we obtain the integral representation for the resolution of the identity operator and, as 
a consequence, for the solution of the Cauchy problem for the master equation as well. 
Then wc show that the expression under the integral can be expanded into a uniformly 
converging power scries, equivalent to the detailed BA, which being integrated term by 
term yields finally the multiple infinite sums coming from the geometric solution and 
from the detailed BA. 

The paper is organized as follows. In section II, we formulate the dynamical rules 
of the model and announce the final result. In section III, we describe the BA for the 
infinite lattice and give the BE for the ring. In section IV, we discuss the details of the 
analytic structure of the BE solutions. In section V we develop a technique of evaluation 
of sums over the Bethe roots. In section VI, we prove a formula for the resolution of 
the identity operator and obtain a formula for transition probabilities in the form of 
infinite sums. Then, to show that the sums are actually finite, we estimate the number 
of nonzero terms they contain. In the last section VII, we discuss the results obtained 
and subsequent perspectives. 

2. Formulation of the model and results 

We consider P particles on the ID ring consisting of L sites. The notation 

P = P/L (1) 

is used for the density of particles. The model has totally asymmetric dynamics, i.e. 
all particles jump only in one direction, which we refer to as forward. At each step of 
discrete time t, a particle from each occupied site. 
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i takes one step forward with probability v 

or 

ii stays with probability {1 — v), 

provided that the target site is vacant. When the next site is occupied, 

iii a particle stays with probability 1. 

All sites are updated simultaneously. We define a configuration X of particles by 
the set of their coordinates X — {xi, X2, ■ ■ ■ ,xp}, written in strictly increasing order, 

Xi < X2 < • • ■ < Xp. (2) 

The finite ring geometry implies also that the coordinates are confined to the values 
from 1 to L. 

^1 > i,xp < L (3) 

The probability Pt{X) for the system to be in a configuration X at time t obeys the 
Markov equation 

Pt+i{X) = J2'^{X,X')P,{X% (4) 

{X'} 

where T{X, X') is the probability of the transition from X' to X for one time step. The 
transition probabihty T{X,X') defined by the above dynamical rules is a product of 
factors, each corresponding to a particular cluster of particles in the initial configuration 
X'. The word "cluster" refers to a group of particles, which has no empty sites between 
occupied sites and two empty sites at the ends. The value of these factors is either v or 
(1 — v) depending on whether or not the first particle of a given cluster jumps during 
the transition from X' to X. 

AfciX') 

T{X,X')= II {1 - v)'-"'' v""' (5) 

i=l 

Here mj = 0, 1 is the number of particles hopping from i-th cluster of X', and jVc i^') 
is the number of clusters in X'. 

It has been shown in [1] that in the case of the infinite lattice, when constraint 
(3) is omitted, the conditional probability, P{X;t\X^;0), for the system to be in 
a configuration X = {xi,X2, ■ ■ ■ ,xp} at time t, given it was in a configuration 
X^ = {x? , Xg, . . . , Xp} at time 0, is the following quotient 

The function JF^^ {X, Y, t) depending on two particle configurations X = {xi,X2, . . . ,Xp} 
and Y = {yi, y2, ■ ■ ■ , yp} and on time t is given by the determinant oi P x P matrix 

■^oo {X, Y, t) = det [f{i - j, Xi - yj, t)]^<ij<p , (7) 
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where the matrix elements are defined in terms of a single function / (a, b, t) expressed 
via the Gauss hypergeometric functions: 



f{a,b,t) = {l-vY< 



( V \ H-t-a)i p I a,-t-a + b 

\v-l) b\ I })+ 1 ' '"^^ 

a — b, —t — a ^ 
-b + 1 



(")-6 jp 



b>0 
b<Q 



(8) 



The notation (a)^ is for the shifted factorial (a)^ = a(a + 1) ■ ■ ■ (a + n — 1). 

The aim of the present article is to show that on the ring of size L the same quantity 
P{X; t\X^; 0) is also given by similar quotient of two terms, both expressed as a single 
function J^l {X, Y, t) of particle configurations X, Y and time t, 

J^L{X,X^,t) 



P(X;i|X";0) = 



(9) 



J'l(X,X,0) 

with the arguments in the numerator and the denominator taken as in (6). However, 
for the finite ring, the function J^l {X, Y, t) is the P-tuple sum of determinants 



J^L{X,Y,t)= J2 ■ 

ni=— oo 



X det 



oo 

np=—(x> 
P 



(10) 



f \ + Prii - 5^ rik, Xi - yj + n^L, t 

' ' l<i,j<P 

unlike the single determinant in the case of infinite lattice. The matrix elements of the 
corresponding matrices are still given in terms of the function f {a,b,t) defined in (8), 
but its arguments depend now not only on the matrix indices i.j but also on summation 
indices Uk- Below we argue that though formally the sums are infinite, they contain only 
a finite number of nonzero terms for any finite time t. Furthermore, the denominators 
of (6) and (9) depend only on the number of clusters in the corresponding configuration 
X being equal to 

xAr,(x)-p 



^oo(X,X,0) = ^^(XX,0) = (l 



Remark 1 One has to define Mc{X) in (11) as a function of X separately for the 
infinite and finite lattices. Indeed, the particles, which occupy the sites 1 and L on the 
ring, belong to the same cluster, while on the infinite lattice they do not. Therefore, 
the value of J\fc{X) can be different for these two cases, even though the coordinates of 
particles formally coincide. 



3. Bethe ansatz. 

3.1. Infinite lattice. 

We first remind the reader of the technique used to deal with the infinite lattice case 
[1]. Consider P particles on the infinite lattice. The particle configurations are given 
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by P— tuples of particle coordinates which are unbounded integers being selected from 
the set 

Z^ = {X eZ^;xi < X2 < ■■■ < xp} . (12) 

Let us introduce the infinite dimensional vector space V^o over the field of complex 
numbers C given by the linear span of the basis 

Xoo = {\X) : X e Z":} , (13) 

i.e. the set of the vectors labelled by the particle configurations. In addition one 
introduces the basis of the dual space V^, which is the span of the dual basis, 

A'^ = {(X|:XGZ^}, (14) 

with the inner product defined by 

{X\X')^S{X,X'). (15) 

Below the bases will be referred to as the configurational left and right bases 

respectively, unlike the left and right eigenbases of the evolution operator to be 
considered. The evolution operator T is defined in terms of the transition probabilities 
T{X, X') defined in (5) 

T= Yl \X)T{X,X'){X'\. (16) 

The problem under consideration is to find the transition probability P(X;t|X°;0) 
from a configuration to X for t steps, which is nothing but the corresponding matrix 
element of the operator T* 

P(X;t|X°;0) = (X|T*|Xo>. 

It was shown in [3], [1] that the evolution operator T has the left and right 
eigenvectors \Bz), {Bz\ parametrized by P-tuple complex parameter 

Z = {zi,...,zp}eC''. 

The eigenvectors corresponding to the same value of Z solve the left and right 
eigenproblems, 

T \Bz) = A{Z) \Bz) , {Bz\ T = A{Z) {Bz\ , (17) 

associated with the same eigenvalue 

p 

A{Z)^{l + X)-^l[{l + Xzi), (18) 

i=l 

where we introduce the parameter A 

The projection {X\Bz) of the right eigenvector \Bz) to the configuration X is given by 
the Bethe ansatz 

(Xi, . . .,Xp\Bz) = W{X)Y,^a,...a,Z;^^' . ..Z-;', (20) 
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supplied with the additional configuration-dependent factor W{X). The latter is 
proportional to the stationary measure of the configuration X and is defined in terms 
of the number of clusters J\fc {X) in the configuration X 

W (X) = (1 + A)^'=^^)-^ . (21) 

The amplitudes A^^, are indexed by the permutations a — {ai, ■ ■ • ,ap} of the integers 
1, . . . ,P. An elementary transposition of two indices (ij) results in the amplitude A^- 
being multiplied by the scattering factor —S{zi,Zj), 

■^...ij... — ~S [zi, Zj) A ji__ , (22) 
of the following form 

^^"-"^■)=1-1/.,1 + A./ ^^^^ 

A remarkable property of the scattering factor S{zi,Zj) is that it is a product of two 
factors, each being dependent only on one of the two parameters Zi, Zj. This property 
allows one to represent the amplitude A^^ in a simple product form 

where each multiple depends only on one of the parameters Zi, . . . , zp. Unlike the right 
eigenvector, the left one, {Bz\ , has no factor Vr(X), while the Bethe part can be 
obtained from that of (20) by the change A^ A~^, Xi —Xi. 

{Bz\x,, ...,xp) = Y, KLaXl • • • C- (25) 

Then, that the projection of the eigenvectors (^Bz\X'^ and {X\Bz) to the basis vectors 
corresponding to a particular configuration X can be represented in the form of 
determinants 

{X\Bz) =VF(X)detB, (26) 
{Bz\X) = detB, (27) 

where the matrix elements Bij and Bij are given by 



1 + \Zi 



^« = ^/^«=(T3I7tj (28) 

One of the main results of [1] is the proof of the formula for the resolution of the identity 
operator. 



/ 



Bz){Bz\dpL{Z)^E, (29) 



where E is the identity operator, dfj,{Z) = {P\)^^ YliLi i.dzi/2nizi) is the integration 
measure, and the integration is performed independently over each Zi along the contour 
r°° encircling the points z = and z = 1, while the point z = — 1/A stays outside 
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Figure 1. Shape of the contour r°°. It must encircle the points z = 1 and z — while 
the point z = — 1/ A stays outside. 

(figure 1). Practically the proof was given by the direct evaluation of the integral in the 
configurational basis, 

j {X'\Bz) {Bz\X) dfi{Z) = 6x,x'- (30) 

As soon as the equality (29) is established, one can insert the identity operator into the 
matrix elements, 

(X|T*|X°> = J {X\T'\Bz){Bz\X'')dfi{Z) (31) 

= y A* (Z) {X\Bz) {Bz\X') dfi{Z) (32) 

so that the final result (6-8) immediately follows from the explicit form of the eigenvalues 
and the eigenvectors. 

One of the observations made in [1] was that there is no any obvious procedure of 
choosing a contour and a measure of integration. However, validity of a particular choice 
can be verified a posteriori by proving the resolution of the identity operator, which is 
fulfilled by a direct evaluation of l.h.s. of (30). It will be clear below that the form of 
the contour can be validated directly by considering the infinite lattice as a limiting case 
of the ring of the size L, which grows to infinity. Then, the Z— spectrum of the model 
is obtained as a continuous limit of the discrete spectrum for the finite system. 
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3.2. The ring. 

What does change when one confines the system to the ring of finite length L? 
Apparently, the total number of particle configurations becomes finite. The same must 
be true for the dimension of vector space Vl defined as the span of the basis 

XL^{\X):XeZ^J, (33) 
where Z^^^ is the domain for particle coordinates on the ring. 

Z^j^ = {X e < xi < X2 < ■ ■ ■ < xp < L} . (34) 

As before, the dual space is spanned by the dual basis 

Xl^{{X\:XeZl^}, (35) 
with the inner product defined by 

{X\X')^5(X,X'). 

Obviously, the dimensions of the spaces are equal to the total number of the basis 
vectors, which is equal to the number of particle configurations: 




dim(T4)=dim(y;)= I . (36) 



The solution of the cigcnproblcm turns out to be analogous to the infinite lattice case. 
The only though important difference is that, due to the finite dimension of the vector 
space, a finite number of independent eigenvectors can exist. Technically this follows 
from the fact that vectors \Bz) and {Bz\ are now the eigenvectors of T only for a finite 
discrete set of the values of the parameter Z. This set is to be defined from the system 
of algebraic equations, which follow from imposing the periodic boundary conditions. 
This is the set one has to sum over, when constructing the resolution of the identity 
operator similar to (29). Of course the latter is correct provided that this set is large 
enough to ensure that the corresponding eigenvectors to form the complete bases of Vl 
and VI 

For further convenience wc slightly generalize the problem. Consider the generating 
function F7(X; t|X°; 0) of the joint probability Pt(X, J; t|X°, 0; 0) for the system to be 
in a configuration X at time t, the total distance travelled by particles being J, given 
the initial configuration X^. 

oo 

F,\X;t\X';0) = {exp{^y,))^ = Y.e''' Pt{X, J; t\X' ,0;0), (37) 

j=o 

The evolution equation for F^' {X,t\X'^ ,0) is similar to the original equation (4) for the 
probability, with the only minor change: the transition probabilities must be multiplied 
by the factor e'"' per each jumping particle: 

mX,X')^ Yl (l-v)^~'^'(e''vr\ (38) 
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Apparently, the limit 7 — > restores the original Markov equation for the probability 
of a configuration 

Pt(X; 0) = lim F^(X; 0). (39) 

7— >o 

All elements of the above Bethe ansatz technique can be directly extended to the case 
of nonzero 7 yielding a minor change in the eigenvalue (18) 

p 

A^ (Z) = (1 + A)-^ J] (1 + e-'Xzi) (40) 



i=l 



and in the scattering factor (23) 



which in turn changes the BA amplitudes, (24) 

^2....,H-i)-n(^)""°'. (42) 

These amplitudes being substituted to the the BA for the right and left 
eigenvectors(20,25) yield an expression for the eigenvectors of T^. All the expressions 
obtained for the infinite lattice are valid for the ring geometry, until the coordinates 
of particles take values at the boundary of the domain (34). For the expressions to be 
valid at the boundary as well, one has to impose periodic boundary conditions 

{xi, X2, . . . , xp} = {^2, . . . , xp, xi + L} . (43) 

Applied to the Bethe ansatz (20,25), the boundary conditions yield a system of algebraic 
Bethe equations (BE) 

^.- = (-1)"-' " 

which fix the spectrum of the parameters Zi^...^zp (for details see [3]). The set of 
solutions Z of (44) defines the sets of left and right eigenvectors 

B={\Bl),ZeZ}, (45) 
{{Wz\,Z^Z]. (46) 

A specific feature of the integrable models is that the same eigenvectors diagonahze 
also a complete set of mutually commuting operators that are as many as the degrees of 
freedom. The simplest example is the translation operator r which translates a particle 
configuration one step forward acting to the right 

T . . . , Xp) = |xi + 1, . . . , + 1) , (47) 

while its adjoint action to the vector of the dual space is the one step backward 
translation 

(Xi, . . . , XpI T = (Xi - 1, . . . , Xp - 1| . 



^ ( TT ^^^^^ ~ ^) + ^^^^0 (AA\ 

^ ' AA^, (^^e^-l) (l + en^,)' ^ > 
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Apparently the vectors from the sets B and B are the eigenvectors of r, 

T\Bl)^Tz\Bl), (48) 
{Bl\r^rz{Bl\, (49) 

with the eigenvalue 

Tz = izi---zp)-\ (50) 
The translation by L steps returns the system to itself, i.e. 

= E, (51) 
where E is the identity operator. As a result tz must be an L-th root of unity. 

(zi---zp)^ = l. (52) 

The same result can be obtained by multiplying all L BE (44). 
Assume that the sets B and B arc dual to each other, 

{Bl\Bl) = {Bl\Bl)5z,z', (Bl|Bl>^0 foranyZ,Z'eZ (53) 

and complete (i.e. their cardinalities are as big as the dimension of the original space 
(36)). Then, the resolution of the identity relation holds 

V ^' = E, 54 

z^z \^z\^z) 

which in the configurational basis reads as follows 

This allows us to derive the matrix element we are looking for. 

Z(.Z \^z\^z/ 

^ (X\B}) (Bl\X^) , , 

- E^7(^) J • (56) 

However the proof of completeness and orthogonality is a separate difficult problem. An 
alternative way, which allows one to obtain the final formula (56) without discussing 
these issues, is to prove the resolution of the identity relation by a direct evaluation of 
the sum on the l.h.s. of (55). In our case, this problem turns out to be solvable without 
explicit knowledge of the spectrum Z. 

4. Location of the solutions of the Bethe equations in complex plane. 

Let us consider the system of BE (44). We introduce for notational convenience a new 
variable 

^nev, ^ ^.g7. (57) 
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Below we will work only with these variables, so we omit the superscript "new" at Zi 
avoiding an abuse of notations. 

Let us gather up into a single constant those parts of the equations (44), which 
have no explicit dependence on the index i, 

--(-r'e-n^^, (58) 

Then the system of BE (44) takes form of a unique polynomial equation of degree L 

z''-'' {z - if - C {1 + zXf = 0. (59) 

For specific values of C, those P of L roots of the polynomial, which match the constraint 
(58), give a solution of the Bethe equations. The formal procedure of finding the solution 
was described in [4], [6]. First, all L roots of the polynomial equation (59) ought to be 
found as functions of the parameter C. Then one substitutes any P of them into the 
equation (58), obtaining a single equation. By solving this equation one obtains the 
solution of the BE corresponding to a given set of the chosen roots. 

Let us consider the analytic structure of the solutions in more detail. One can 
rewrite (59) in the form 

w{zf = C, (60) 

where the function w{z) is 

^i-p (z - lY 

Equivalently one can write 

w{z) = exp {2mk/L)C^/^, (62) 

where the integer k is any integer chosen from the range 1 < k < L specifying a 
particular choice of the branch. The branch of C^^^ is implied to be fixed, e.g. 

< arg(CV^) < ^. 

We use the notation z+, Z- for the two solutions of the equation 

dw (2;) 



dz 

which yields 



= 0, (63) 



-1 + A - 2Ap ± V(l + A)(l + A(l - 2p)^)) 

^ ^ "w^- ■ ^ ^ 

We define the domain D^, (figure 2a), as the extended complex plane cut along 
the segments [0,1] and [— 00, — 1/A] of the real axis and punctured at the points 
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Figure 2. The domains Bz (a) and Du, (b). 



{0, 1, — 1/A, oo, z_}, and also the domain D^, ( figure 2b), as the extended complex 
plane punctured at the points 0, oo, and at four points 

wl = \w{z+)\e'''P, (65) 

wl = \w{z+)\e^'''-''"', (66) 

wl = |«;(^_)|e"'-'"^ (67) 

^2 = \w{z^)\e^'+'''P, (68) 

and cut along the straight segments: [0,w^] , [O,?/;^ , [wi,oo] , [wi,oo]. If we define 
©2 as a domain of w{z), such that the value of arg [^(-z)] is vrp, vr (2 — p), and tt (1 — p) , 
TT (1 + p) at the upper and lower banks of the first and second branch cuts of 
respectively, then the mapping w{z) : — > is the monovalued analytic mapping. 
Going around the branch cut [0, 1] the value of arg [w(z)] changes by 2n exactly once. 
Therefore, no repetition in the value of w{z) can occur, i.e. the mapping w{z) is schlicht 
(one-sheet). Hence the mapping w{z) is bijective and one can construct an analytic, 
monovalued, schlicht mapping u'^^(-) : —>■ D^, inverse of w{z). Then, given a complex 
number a G D^, the equation w{z) = a has a unique simple root z{a) G D^, being an 
analytic function of a in D„,. 

Thus, given a value of C, each root of (60) is the unique simple root of (62) for 
some k, and all the roots are different for different k-s. To solve the Bethe equations 
one formally can choose P integers ki, . . . ,kp from the set 1, . . . , L. Then, substituting 

zj = vj-^ [exp {27rikj/L) C^l^\ (69) 

into the constraint (58) for all j = 1, . . . , P, one obtains a unique equation for the 
parameter C . Solving the equation for C, we obtain the solutions corresponding to 
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given set ki, ... ,kp. Going through all the possible sets of integers ki, . . . ,kp we obtain 
all the solutions of the Bethe equations. 

It is easily seen that the eigenvectors (20,25) identically vanish on the solutions 
corresponding to the sets ki, . . . ,kp containing a pair of two equal integ ers, ki — /tj, 
which result in Zi — zj. Therefore one has to look over only those sets ki, . . . , kp, where 
all the numbers are different. If we assume that for any set the equation for C has 
exactly one solution we obtain just as many solutions as we need (36) to obtain the 
complete set of hnearly independent eigenvectors. The direct proof of this fact however 
is beyond the aims of present article and will be considered elsewhere. 

Remark 2 For generic values of the parameter 7, and hence of the parameter C, the 
r.h.s. of (62) is away from the branch points ofw^^{z), which ensures that the root of 
(62) is simple. However, for specific values of ^ this can he not the case. This does not 
create a problem as this situation can be considered as the limiting case of the generic 
one. In particular, a double root can appear at z+ or z^, as a consequence of square root 
singularities of w~^{z). In this case, one can think of it as a pair of roots at different 
banks of the branch cut. This in fact specifies the way in which it evolves when the 
value of J changes. Also the limit C ^ 0, which implies 7^0, corresponds to P roots 
meeting at z = 1. Then, the choice of different integers ki for different roots Zi removes 
the degeneracy as the arguments of {zi — 1) are different. 

Let us consider the curve Fc defined by the equation 



where c is a real, positive number. All the roots are located on curves where c — |C^''^|, 
for some discrete set of values of C related to the solution via the equation (58). The 
particular case A = of was named in [6] as the generahzed Cassini oval. To describe 
the form of this curve, let us first look at the behaviour of \w {z)\ at the real axis. It 
is easily seen that \w {z)\ has two zeroes at the points 2; = and -2 = 1, and diverges 
for z ±00 and z — 1/A. There are also two extremums: z^ and z_ given by 
(64), which are a minimum and a maximum respectively. As p varies from to 1, the 
point Zj^ monotonously moves from 1 to and z_ from —1/A to —00. At the same 
time |w(2;+)| increases from zero for < p < 1/2, reaches the maximum at p = 1/2 
and then decreases back to zero for 1/2 < p < 1, while (2;_)| decreases from infinity 
for < p < 1/2, reaches the minimum at p — 1/2 and then increases back to infinity 
for 1/2 < p < 1. Note that the values of \w (^+)| and \w {z-)\ meet only in one point 
p — 1/2 approaching its values from below and above respectively. Let us look at the 
cross points of the plots y = \w (x)| and y = c as the constant c grows starting from 
zero. The following stages exist (see figure 3): 

c < 1^(2;+)! For small c there are four cross points {—1/A < < 0, < 2:2 < -^3 < 1, 
1 < Z4}, which corresponds to two "ovals" encircling the origin and the point z = 1. 
In the limit c ^ the form of the contours approaches circles of vanishing radius 




(70) 
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Figure 3. Schematic picture of tlie stages of the evolution of the contour Tc as c grows 
from zero to infinity. At the lower plots the the graphs y = \w{x)\ and y = c are shown 
for each stage, which cross points define the cross points of Fc with the real axis. 

collapsing to the points z = 0,1. As the value of c increases the point Z2 and 
move towards each other, and the radius of the "ovals" increases . 

c = When the value of c reaches the points Z2 and z^ merge at z+, i.e. 

the two "ovals" develop cusps meeting at 2+. Then the shape of the curve resembles 
the lemniscate. This form survives, when one considers the thermodynamic limit. 
Specifically, the right part of the curve is the one considered in studies of low lying 
eigenstates [3]. 

|tf(24.)| < c < As c exceeds |w(2;+)| the crosspoints ^2,^3 disappear after merging 

at 2;_|_such that only the two, { — 1/A < < 0, 1 < z^} remain. These are two 
crosspoints with the horizontal axis of the "big oval", which appears after the 
two "smaller ovals" merge. The big oval contains the points z = 0,1 inside, while 
the point z = —1/A stays outside. As c increases zi,Z4 move along the real axis 
towards —1/A and +00 respectively. As the point zi is confined between and 
— 1/A, while the effective radius of the "oval" is unbounded, it finally starts to bend 
around —1/A and approaches the real axis from below and above at the right of 
the point —1/A. 

c = At this stage the two points at the parts of the "oval" bending around —1/A 

meet at z^, so that two "ovals" appear, one inside the other. They have sharp cusps 
at their only common point z_ . The bigger "oval" goes around —1/A, 0,1 while 
the smaller encircles only —1/A. 
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c > After forming two ovals at previous stage, they detach from each other by 

sphtting the common cross point with the real axis at z_ into two new cross points 
Z2, z'^, which then move along the real axis to —00 and — 1/A respectively as c go 
on growing. The form of the contours then approach two circles of infinite and zero 
radia, the latter collapsing to the point —1/A. 

In the case p = 1/2, when the values of \w {z+)\ and \w {z^)\ are equal, the second 
and the fourth stages coincide while the third one does not take place. 

An important point of the above analysis is the limiting shape of the curves under 
consideration as c goes to zero or infinity. Specifically, one can always choose a value 
of c so small that the contour Fc separates a small neighborhood of the points and 1 
from the rest of the complex plane. For c large, the same is valid for the points —1/A 
and 00. This fact will be used in the next section to evaluate the sum over the Bethe 
roots. 

5. Summation over the Bethe roots. 

Let us rewrite the Bethe equations in terms of the variables {zi}, given by (57), in 
polynomial form 

Pi{Z) = 0, (71) 

where it is convenient to write the polynomials Pi{Z) as a difference of two other 
polynomials 

P,{Z) = gi{Z)-hi{Z) (72) 
which read as follows 

p 

g, (Z) = e-^^^^^ {z, - if l[z,{l + z,X) , (73) 

i=i 

K (Z) = (-1)^-^ (1 + z,Xf f[ (zj - 1) . (74) 

The aim of the present section is to evaluate the sum of the analytic functions over the 
roots of the system (71). Two remarks are necessary. First, we imply that the roots are 
bounded from infinity, as otherwise the eigenvectors and eigenvalues would be singular. 
Second, the polynomial equations can be satisfied with the solutions constructed from 
the roots from the set {0, —1/A, 1}. Furthermore, it is easy to see that if one of the 
roots Zi is taken from this set, the other P — 1 roots must belong to this set as well. 
However, solutions constructed in this way obviously do not match the constraint (52) 
for 7 7^ 0, and, therefore, must be excluded. Appearance of extra solutions is due to 
multiplication of the BE by an expression that itself can be zero, when transforming 
it to the polynomial form. In the case 7 = one such a solution exists. Namely it is 
Zi = ■ ■ ■ = zp = 1, which corresponds to the ground state of the evolution operator, 
i.e. the stationary state of the stochastic process. As will be seen below, this case can 
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be treated as a limiting case of the generic situation and does not require a special 
consideration. 

Let us define the domain D C C of the complex plane by the inequalities 

D^{z :a<\w{z)\<A}, (75) 

where a and A are two real positive constants, such that a < A. Denote by Fc the 
contour discussed in the previous section, defined by (70) 

r, = {z:\w{z)\=c}. (76) 

Then the domain D is between Ta and Ta, and its boundary 

dD = {z:TaUTA} (77) 

is oriented in such a way that going along it in positive direction one keeps the interior 
of D left. Then we form a polycylinder domain Q cC^ as a cartesian product of P 
copies of D 

DiX ■■■ X Dp. (78) 

The skeleton F of Q is the subset of its boundary dfl consisting of points, which are at 
the boundary of every Di, . . . , Dp: 

F = 9Di X • • • X dDp. 

The definition (75) of D guarantees that all the points Z = {zi, . . . , zp), such that 
Zi e {0,— l/A,l,oo} for some i, are outside of Q. On the other hand, the dimension of 
the complement of Q approaches P — 1 as a and A go to zero and infinity respectively, 
i.e. their Lebesgue measure in vanishes. Thus, it is natural to expect that for a 
small and A large enough all the roots of the system (71) fall into Q. Then, the sum of 
an analytic in function f{Z) over the roots of BE which fall into Q can be evaluated 
with the aid of the multi-dimensional logarithmic residue theorem [9] . 

Theorem 3 Let Qc be a polycylinder domain with piecewise smooth boundary dfl 
and F be its skeleton. Let the mapping {Pi{Z),i — 1, . . . , P} : ft — > be holomorphic 

in O and have no zeroes at the boundary dfl. Then, for any function f{Z) analytic in 
ft, the sum of its values over the set Z — {Z & ft -.Pi {Z) — 0,i — 1, . . . , P} is given by 
the following integral 



t(7\ 1 / ..^. dP.iZ) , , dPp(Z) 



det 



r nf=i Pi {Z) 



dPj jZ) 
dzj 



p 



l<i,3<P 1=1 



Every zero is counted as many times as its multiplicity is. 

The theorem is a particular case of the CaccippoUy, Martinelly, Bishop, Sorani 
theorem, see [9]. The original theorem is proved for the domain called special analytic 
polyhedra, which in particular ensures the absence of zeroes of the mapping at the 
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boundary of domain, which in turn guarantees that all zeroes inside the domain are the 
isolated ones. In our case the absence of zeroes at the boundary dfl of the domain Q, 
is provided by the following lemma. 

Lemma 4 Let fl be defined as above. Let A < 1. Then in the range of 7 : 
(p~^ |liiA| > 7 > 0), there exist constants Qq and Aq, such that for any a < Qq and 
A > Aq the mapping {Pi{Z), . . . , Pp{Z)} has no zeroes on dfl. 

Proof. Suppose there is a point Z — {zi,...,zp} e dCl, such that Pi{Z) — for 
i — 1, . . . , P. It implies that 

9^{Z) - ----- 



K{z) 



= e 



-1,7 



wyzi 



Zi---Zp 



1. 



(80) 



w{zi) ■ ■ ■w{zp) 

The point Z being at the boundary 917 of the domain 17 means that at least one 
coordinate Zi from the set ^i, . . . , 2;p is at the boundary of D^, i.e. is either on or F^. 
Let first 



Zi e Fa 



(81) 



for some i. As Zi enters into (80) only via \w{zi)\^, while the other factors do not 
depend on the index i at all, \w{zi)\^ does not depend on i either, i.e. (81) holds for all 



1, . . . , P, which immediately yields 



Zi--- Zp 



12) 



Recall that a always can be chosen small enough such that Fq belongs to small 
neighborhoods of the points z = Q and z = 1. In other words for any small £ > 
one can choose a small Oq such that for any a < Qq 

(83) 



sup 



^ < 1 + e. 



Taking e < {e'^^ — 1) we obtain 

ki ■ ■ ■ Zp 



^1" < e^\ 



(84) 



which contradicts (82). 

Consider now the case when Zi G F^ for some i. Analogously to the previous case 
one needs to satisfy (82) with the set of solutions 2:1, ... , 2;p, which are on F^, i.e. either 
go to infinity or to (—1/ A), as A increases. Thus, for any £ > one can always choose 
a large Aq such that for any A> Aq 

1 



inf 

{z&Ta} 



z\ > 



A 



— e. 



(85) 



Taking £ < (|1/A| - e^P) we obtain 

\zi---zp\^'p > e^"-, (86) 

which also contradicts with 82. ■ 

This lemma in particular ensures that the constants a and A always can be chosen 
small and large respectively, so that all the solutions of BE, are inside 17. Indeed, 



Determinant solution for the TASEP with parallel update II. Ring geometry 



19 



according to the previous section the solutions always belong some contour Tc for a 
particular value of c. Hence they would belong to the boundary of f2 defined with c — a 
OT c = A. As follows from the above lemma no such solutions exist for a < oq and 
A > Aq. Therefore all the solutions are in f2 defined with such a and A. 

The next lemma provides a condition for evaluating the integral in (79) in the form 
of infinite series. 

Lemma 5 Let the domain O and its skeleton T be defined as above. Let the range of 
parameters be 



and 



< A < 1, 

|ln A| < 7 < 

1 

o< -. 



(87) 
(88) 

(89) 



Then, one can choose the constants a and A such that for any Z E T the following 
conditions hold: 

9i{Z) 



if Zi e Fa and 



hi{Z) 



9i{Z) 



hi{Z) 



< 1 



> 1 



(90) 



(91) 



if Zi e Furthermore the limits a — > and A — > oo can be taken simultaneously, in 
such a way, that these inequalities hold. 

Proof. It follows from the explicit form of w{z) that for any £ > one can choose a so 
small that 



inf \z\ > (1 — e) a^-p 
2er„ 

sup \z\ < {1 + e) . 
On the other hand for any £ > one can choose A so large that 

(1-6) 



inf l^;! > 



A 



(92) 
(93) 

(94) 
(95) 

(96) 

Now we can use the inequalities (92) and (94) to estimate the bounds for the ratios 
\zi/ci\, which yield 

< ^^"'^P < (l + £)max((AA)T^,^ ) .(97) 



sup 1^1 < (1 +£) Ai-pAi-p. 

Let us write Ci = a if Zi E and Ci = A if Zi E Ta- Then we have 



g^iz) 




Zi - ■ ■ Zp 






Ci ■ ■ ■ Cp 



;i -£)min ( 



Ci • • • Cp 



Determinant solution for the TASEP with parallel update II. Ring geometry 



20 



If we choose A and a such that 



then for c,- = a we have 



-L-y L 



Zi---Zp 



<(l + £)^e-^^<l, 



Ci ■ ■ ■ Cp 

which is (90). At the same time if we choose them such that 
^ 1 

then for c,- — A we have 



2;i • ■ • 2;p 



Ci • • • Cp 



1/p 



> e 



-L7 



1 -e 
A 



>1, 



(98) 



(99) 



(100) 



(101) 



which is (91). The only what we need is to satisfy both conditions (98) and (100) 
simultaneously, which implies 

(XA)-^ < a w < (AA)-(^)' . (102) 
This is possible if p < 1/2. ■ 

Remcirk 6 The limit 7 ^ m two above lemmas can be considered after the limits 
a — > and ^4 ^ 00 are taken. This particularly solves the problem of the groundstate 
mentioned above. While for 7 > the roots corresponding to the groundstate are away 
from Z = (1,...,1), they approach this point when 7 approaches zero. However due to 
the order of limits described they still remain separated from this point by the boundary 

ofn. 

Using lemma 5, we can represent the factor l/Pi[Z) under the integral in (79) in 
the form of infinite sum. Indeed, for Zi e we have 



1 



1 



P^{Z) 

while for Zi e Fa 

1 



gi (Z) - K {Z) hi (Z) 



n=0 



Pi{Z) gi (Z) - hi (Z) hi (Z) 



E 



9dZ)_ 
hi{Z) 



hi{Z)) 



(103) 



(104) 



Thanks to (99,101), the series are absolutely and uniformly convergent, and, as such, 
can be integrated term by term. Note that the summands have no singularities in 
f2. Therefore the contours and F^ can be deformed into a single contour. With 
respect to singularities of the expression under the integral, it has the same form as the 
contour F°° described in the Section 3, which was used to construct the resolution of 
the identity operator in the case of infinite lattice. Thus, both the contours F^ and Va 
can be deformed into F°° and the term corresponding to the integral over F^ brings a 
minus sign due to the opposite orientation. The final expression has no dependence on 
the values of a and A. Therefore, assuming the limit a — > 0, A — > cxo one can think of 
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fl as of product of P complex plains punctured in four points {0, 1, — 1/A, oo}, which 
clearly contains all the necessary solutions of the BE. As a result we have the following 
expression for the sum over the roots of the Bethe equations. 

Theorem 7 The sum of the values of a function f{Z) analytic in the whole complex 
plane, except maybe the points {0, 1, —1/A, oo}, over the roots of the Bethe equations is 
given by the following P-tuple absolutely converging sum. 

'hi{zy 

/ poo v ... V T"^ 

np=-oo^ ^ 1 ^ P 

"cl(lngi(Z)-lnWZ))' 



E/(2)= E ■•• E 



9i{Z) 



X det 



n 



J l<k,j<P 1=1 



dzi 
2m 



(105) 



Note that we write d {In gk{Z)) /dzj = {dgk{Z)/dzj) / gk{Z) in the determinant 
instead of {dgk{Z)/dzj) /hk{Z), because the equality hk{Z) = gk{Z) holds on the roots 
of BE, which are the only contributing the integral. 



6. Proof of the resolution of the identity and formula for the transition 
probability 

Now, we are in a position to write the sum in (55) in the integral form (105) substituting 

To this end we need the expression for the norm (B^ze-yl^ze—y)- "^^^ hypothesis about 
the form of the norms of Bethe vectors was first proposed by Gaudin [10, 11]. Later 
it was proved by Korepin [12] within the quantum inverse scattering method for XXX 
and XXZ type models. Those results can be applied to our model with minor changes. 
However, they require developing the quantum inverse scattering method, which is not 
a subject of the present article, and will be done elsewhere. Here we use the formulas as 
they are given in the Korepin's article. Note that the proof below and of the final result 
do not rely on the validity of the formula for the norm. The latter serves as a hint for 
writing the expression under the integral, while the validity of the results follows from 
the resolution of the identity proof made independently. Written in the the transformed 
variables (57), Gaudin formula yields 



(^e-.l^le-.>=det 



(107) 

i,j=i,...,p 



' dzi hj (Z) 

which being in the denominator of the expression under the integral cancels the Jacobian 
in the numerator, yielding only the factor {zi • • • zp)~^. Substituting the explicit form 
of the eigenvectors (20,25) and the functions g{z) and h{z), (73,74), and performing one 
summation over the permutations, which is trivial due to the permutation symmetry of 
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the summands, we come to the following expression. 



Z^Z \ ^' 2 1 ni=— oo np=— CX) <tGSp 



X 



/ 27ri2;i 
Though there is an infinite P-tuple sum, it turns out that only few terms in each sum 
contribute. The following lemma establishes which summands are nonzero. 

Lemma 8 Let X,Y e be two particle configurations, a be a permutation 

(ui, . . . , CTp) of the integers 1, . . . , P, {ni, . . . ,np} e be a set of integers. Then, 
the necessary conditions for the product 

to be nonzero are 

n^G {-1,0,1} (110) 

for i = 1, . . . , P. Furthermore, the cases when Ui = ±1 for some i = 1, . . . , P, suggest 
that 

p 

1=1 

and there are clusters in the configurations X and Y, which in particular contain the 
sites 1 and L. 

Proof. To evaluate the integrals over under the product we expand the expression 
under the integral into the Laurent scries in the ring 1 < |2;| < 1/A, implying A > 1, 
and look for the coefficient coming with z^^. Let us for brevity introduce the notations 

p 

dij ^ i - j + Pui-^ rife, (111) 
k=i 

bij = Xi- yj + Lui. (112) 
Then, the integrals to be nonzero, the following conditions must be met: 

a. if Ojo-. = 0, then = 0; (113) 

b. if ajo-. > 0, then 6^^. < Oj^.; (114) 

c. if ttiai < 0, then > a^^.. (115) 

Given particle configurations X, Y and a permutation o", every element from the set 
{?T'i}i=i,...,p fall into one of three classes (a),(b) or (c) depending on the sign of Ojo-.. 
Then the inequalities (113-115) determine which restrictions on Ui must be imposed, all 
the integrals in the product to be nonzero simultaneously. 

First we note that as the coordinates of particles on the finite lattice satisfy 
\xi — yj\ < L, the only way to satisfy the equality bi^- = is to put Ui = 0. Therefore, 
this is always the case for the numbers n^, which belong to the class (a). 
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One can use a similar argument to consider a particular case of the set {ni}i=i^,,,^p 
with all components equal. 

rii — 712 — ■■■ — np = n. (116) 

Then for such i that Ojo-i — {'i — (^i) > 0, class (b), (114) requires n < 0, while for such j 
that ajfj. — {j — (7j) < 0, class (c), (115) implies n > 0. Apparently for any permutation 
(T, where for some i there is a number cTj which satisfies (i — Ui) > 0, there must be at 
least one j such that {j — (Tj) < 0. Thus we necessarily have n — 0. 

To study the other sets we estimate the bounds of their maximal and minimal 
elements. Consider the sets {^i}i=i,...,p, where not all components are equal. From such 
set a maximal element can be chosen 

JT-max = max{ni, . . . , Up}. (117) 

Of course, in general several numbers from the set can attend the maximum. The proof 
below consists of two steps. We first show that at least one of them falls either into 
the class (a) or into the class (b), i.e. satisfies either the first equality in (113) or the 
first inequality in (114). Then we find which restrictions on nmaxfoUow from the second 
equality of (113) or the second inequality (114). 

To proceed with the first step we note that for any k such that = rimax at least 
a part of akak is always positive 

p 

Pnfc-^ni>0. (118) 

1=1 

The remaining term {k — ak), which also enters into afeo-^,, is either positive or negative. 
In the case when it is positive or it is negative, but its absolute value is smaller then one 
of the other terms, we have Ofco-j, > 0. It can happen, however, that {k — ak) is negative 
and its absolute value is large enough to turn aka^ to be negative. This is possible when 

p 

Pnk-J2nk<\k-ak\<P-l. (119) 

k=l 

Let s be an integer, 1 < s < P, such that 

p 

Pnu.^^-Y,^k^P-s. (120) 

k=l 

The last equation suggests that the number elements in the set {Tii}i=i,...,p equal to nmax 
is at least s, i.e. there exists at least s integers 1 < /ci, < • • • < /c^, < P such that 

J^fel = • • • = JT-fe. = JT'max- (121) 

Let us now check the sign of aic^ for i chosen from ki, . . . ,ks and an arbitrary permutation 
a. At worst we have = P — s + 1, . . . , (Xfc^ = P, in which case for i — ki, 
Qifj^ = [ki — {P — s + 1)] + (P — s) > 0. For any other a either cTj can be taken 
smaller or i bigger, so the inequality becomes strict. Thus, we conclude that there is 
always at least one maximal element belonging either to the class (a) or (b), which 
finishes the first step of the proof. 
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According to the above arguments for the class (a) we get Umax = 0. Let now 
rii = nmax be in the class (b). Combination of the second inequality (114) and the finite 
lattice restriction on the difference of the particle coordinates 

\(x,-yj)-(i-j)\<L-P (122) 

give the upper bound for ^max- 

rimax < 1 - ^ _ p ^ nk. (123) 

k=l 

All the above arguments can be applied also to the minimal element as well 

rimin = min{ni, . . . ,np}, 

so we conclude that it belongs either to the class (a), where rimin = or to the class (c), 
where (115) and (122) yield 

nmin > -1 - rik. (124) 

fc=l 

It follows from (123) that if ^f^^rifc is positive, then rimax must be strictly less than 
1 and hence nonsensitive. At the same time, one sees from (123) that if X^f^i^ifc is 
negative then nmin must be strictly greater than —1 and hence nonnegative. Thus the 
only possibility is to have 

p 

Y^nk^O. (125) 

k=l 

Substituting this back into (123,124) we conclude that there are only three possible 
values of the elements of the set ,,,^p 

{-1,0,1}, i = l,...,P, (126) 

the numbers +1 and —1 appearing in the set {nj}j=i,...,p equally many times. 
Consider now the case when not all rii are equal to zero, i.e. 

rimin = -1, rimax = 1 (127) 

If we return to the inequahties (114,115) and repeat the derivation of (123,124) using 
(125), we obtain that the case (127) is reahzed when the weak inequality (122) turns 
to the equality. Specifically, deriving the estimate for rii — n^ax — 1, the necessary 
condition for ^max to be equal to 1 is 

Vai -Xi-(ai-i)^L-P. (128) 

This means that the particles with coordinates Xi and belong to the clusters in X 
and Y, which spread to the last and the first sites respectively, i.e. the particle at Xi 
belongs to the cluster of X, which starts with xi — 1 and the particle at y^, belongs to 
the cluster of Y, which ends with y^ — L. Similarly, for rij — rimin = — 1 we have 

- Xj - {aj -j)^-{L-P), 



Determinant solution for the TASEP with parallel update II. Ring geometry 



25 



i.e. the particle at Xj belongs to the cluster of X, which ends with — L and the 
particle at i/o-j belongs to the cluster of Y, which starts with yi = 1. Thus wc conclude 
that in both X and Y there exist the clusters, which cover at least the sites at the 
positions 1 and L. This proves the last statement of the lemma. ■ 

Of course the choice of the reference point at the ring has a conventional character. 
One can get reed of the terms rii — ±1 by simple rotation, which places any hole of one 
of the configurations X or Y either into the site 1 or L. This fact is useful for the proof 
of the resolution of the identity relation. 

Theorem 9 The resolution of the identity operator is given by l.h.s of (54)- 

Proof. Consider the integral representation (108) of the sum in (54). According to the 
pervious lemma if one of the configurations X and Y are such that any of the sites 1 and 
L is empty the only term of the sum that remains corresponds to ni = ■■■ = np = 0. 
This term coincides with the resolution of the identity operator for the infinite lattice. 
For the proof of the infinite lattice case we refer the reader to the Proposition 1 from 
our first paper [1]. 

The cases where in both configurations X and Y there is a cluster containing the 
sites 1 and L can be reduced to the previous situation by translation. Specifically we 
note that the product {X\Bl) (!B^|F) is invariant under the translations, i.e. 



Hence we can repeatedly apply the translation operators r and until a hole comes 
either to the site 1 or L. In this way the problem is reduced to the proved one. ■ 

Since wc have proven the formula for the resolution of the identity operators, we 
can apply it to find the matrix clement we are looking for. To this end, as shown in (56) 
we must insert the eigenvalue of T* under the integral. Then the integral representation 



{xi + l,...,xp + l\Bl) (Bllyi + l,...,yp + l) 
{Xr-'\Bl) (bUtY) = TzT^' {X\Bl) {bI\Y) 
{xi, xp\Bl) (bIIvi, ...,yp). 



of F7(X,t|F,0) is 

(X|T^y> = (1 + A)-^*PF(X) ^ (-1)''^' 



cr£Sp 




i — cji+t 



hi {Z) dzi 



Pi {Z) 2niZi 



(129) 



oo 



oo 



{X\Tt^Y)^W{X) E E ^' 




X 



(130) 



where 




27ri^ 



dz 



(131) 
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The integral for /(a, b, t) is evaluated in terms of the hypergeometric functions (8). The 
sum over the permutations leads us to the determinant, given in (10). By putting 7 = 
we obtain the result for conditional probability announced. 

Like the sum obtained for the resolution of the identity operator, which is the 
particular case of the sum (130) at t = 0, the latter sum being formally infinite, however 
contains finitely many nonzero terms. The analysis similar to one of the lemma 8 shows 
that at time t the upper bound for the maximal element ninax(^) of the set {nj}j=i p, 
which ensures corresponding summand to be nonzero, is 

1 ^ t 
nr.^{t) < 1 - ^^3^ + (^3^' (132) 

while that for minimal one, nniin(^), is still like it was for the t = case 

rirain {t) > -1 - _ ^ Wfe. (133) 

k=l 

The bounds for the sum ^^^^ Uk can be obtained from the following arguments. Suppose 
that X^feLi'^fc ^- Then (132) requires nniax(^) < 0, which contradicts the assumption. 
Thus we have 

p 

J]nfc<t (134) 
k=i 

Another argument can be given, based on the fact that ^k=i^k < -P'T'max(^)- Then 
using (132) we obtain 

fc=i 

The first upper bound is lower than the second, when t < P, and vice versa when t > P. 
For the lower bound we suppose that ^^^=1 nk < 0, which contradicts (133) and results 
in 

p 

J]nfc>0. (136) 

k=l 

Finally we can use these inequalities to estimate the range of the summation indices 
corresponding to the summands, which give nonzero contribution. 

Lemma 10 For the summands of the sum (130) to he nonzero it is necessary that 

""'-^^ {L-py ^"""-^ ^^^^^ 

and 

f P 
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if P <t and 



ni> -1- 



P{L-P + t) 
L{L-P) 



r. t-P 
, Y,nk<P+^- 

k=l 



(139) 



if P > t, for i = 1, . . . , P. In the cases, when the expressions in r.h.s. are not integer, 
the inequalities are strict. 

From the definition of tlie generating function Ff{X, t\Y, 0) one concludes that the 
coefficient of e'^'^ for some nonnegative integer J is the probability Pt{X, J;t\Y,0;0) 
for the total distance travelled by particles for time t and the final configurations X, 
given the initial configuration is Y. One can see that in (130) the similar term is 
Q'Ei=i{^i-yi+Lni) ^ while the probability for the travelled distance to be J is the sum of 
the coefficients of terms, where the sum ni + ••• + np is fixed, X^^^ {xi — yi + Lni) — J. 
Thus, the sum (ni + • • • + np) has a meaning of the total number of windings around the 
lattice all the particles made. By this reason, the sum is always nonnegative unlike the 
individual numbers ni, . . . ,np. The meaning of the latter is well understood in frame of 
the geometric approach to the BA [1]. These are the winding numbers of "virtual" free 
trajectories, which being weighted with corresponding weights can be used to reconstruct 
the TASEP dynamics. 

7. Conclusion and discussion 

To conclude we have obtained the probability of the transition from one configuration 
to another for arbitrary time for the TASEP with parallel update on a ring. To this end 
we developed the method of summation over the solutions of the Bethe equations, which 
is based on the multidimensional version of Cauchy residue theorem. In this way the 
integral representation of the solution is obtained. The expressions under the integral 
can be expanded into the uniformly convergent power series, which being integrated 
term by term, yields the result in form of multiple and formally infinite sum of the 
terms, each having the determinant form. It is shown that only finitely many terms of 
this sum are nonzero. Note that though the convergence of the series under the integral 
is proved for the domain < p < 1/2, the behaviour of the final finite sums have no 
singularities at the point p = 1/2 as well as for any values of p. Therefore we expect 
that arguments of analytic continuation exist which extend the proof for any value of 
density, < p < 1. On the other hand the case l/2<p<lis related to0<p<l/2by 
particle hole symmetry. It is an interesting exercise to find an explicit relation between 
the final formulae of the transition probabilities for these two cases. 

There are several directions of possible development of the result. First, it looks 
possible to generalize the method to extract not only the sum over the solutions of 
the BE, but also to extract the contribution of particular solutions. In this way 
one could obtain closed exact expressions for particular eigenvalues and eigenvectors, 
rather than only the asymptotic behaviour studied before. Second, the integral 
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representation obtained can be useful to study the large time asymptotics for the growth 
phenomena with time. Many similar results where obtained recently for the infinite 
lattice, due to the observed parallels with the theory of random matrix ensembles 
[13, 14, 15, 16, 17, 18, 19]. The Bethe ansatz, giving the integral representations of 
the physical quantities like particle current probability distribution, could also be a 
starting point of such an asymptotical analysis. Particularly, the result of present paper 
could be used to make an advance for the ring geometry where not much results have 
been obtained yet. 
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